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Abstract 

The common view ttiat structure functions measured in deep inelastic 
lepton scattering are determined by the probability of finding quarks and 
gluons in the target is not correct in gauge theory. We show that gluon 
exchange between the fast, outgoing partons and target spectators, which 
is usually assumed to be an irrelevant gauge artifact, affects the leading 
twist structure functions in a profound way. This observation removes the 
apparent contradiction between the projectile (eikonal) and target (par- 
ton model) views of diffractive and small xb phenomena. The diffractive 
scattering of the fast outgoing quarks on spectators in the target causes 
shadowing in the DIS cross section. Thus the depletion of the nuclear 
structure functions is not intrinsic to the wave function of the nucleus, but 
is a coherent effect arising from the destructive interference of diffractive 
channels induced by final state interactions. This is consistent with the 
Glauber-Gribov interpretation of shadowing as a rescattering effect. 
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1. Introduction 



Deep inelastic lepton scattering, C.N — ^ r + X (DIS) is central for our understand- 
ing of hadron structure. Ever since the earliest days of the parton model, it has 
been assumed that the leading-twist structure functions Q^) measured in 
deep inelastic lepton scattering are determined by the probability to find quarks 
and gluons in the target . This probability is given by the target wave function 
at the light-cone (LC) time when the current interacts (in the < frame). 
For example, the quark probability is distribution is 



Y[ dxi (fk^ 



\i^n{Xi,k±i)\'^^S{xB - Xj) (1) 



J=<] 



where the ipn are LC wave functions of the target (see Eq. (Q) below). The 
identification of structure functions with the square of light-front wave functions 
is usually made in the ghost-free LC gauge n ■ A = A'^ = 0, the argument being 
that the path-ordered exponential in the operator product appearing in parton 
distributions (see Eq. (^) below) reduces to unity. Thus the DIS cross section 
appears to be fully determined by the probability distribution of partons in the 
target. 

However, we shall show that this parton model interpretation of the structure 
functions, which was established for a theory with Yukawa couplings [0], is not 
correct in gauge theory. The critical issue is whether the scattering taking place 
after the virtual photon interacts can affect the leading twist cross section. It 
is well-known that in Feynman and other covariant gauges one has to include 
corrections to the "handbag" diagram due to final state interactions of the struck 
quark with the gauge field of the target. Light-cone gauge is singular - in par- 
ticular, the gluon propagator dl^^k) = -g^^" + "^^"^'^fe''"'' has a pole at 
fc+ = which requires an analytic prescription. In final-state scattering involv- 
ing on-shell intermediate states, the exchanged momentum k~^ is of (9(l/z/) in 
the target rest frame, which enhances the second term in the propagator. This 
enhancement allows rescattering to contribute at leading twist even in light-cone 
gauge. 

We find that gluon exchange between the outgoing quarks and target spectators, 
which is usually assumed to be suppressed in the Bjorken limit, affects the leading 
twist structure functions in a profound way. Final state diffractive scattering gives 
rise to interference effects in the DIS cross section. Thus nuclear shadowing is 
not caused by the wave function of the nucleus, but is induced by final state 
interactions. 

Thus the depletion of the nuclear structure functions is not intrinsic to the wave 
function of the nucleus, but in fact is a coherent effect reflecting the destructive 
interference of diffractive channels induced by the flnal state interactions. 
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The distinction between structure functions and parton probabilities is already 
implied by the Glauber-Gribov picture of nuclear shadowing 0, H, |, ^. In this 
framework shadowing arises from interference between complex rescattering am- 
plitudes involving on-shell intermediate states. In contrast, the wave function of 
a stable target is strictly real since it does not have on energy-shell configurations. 
A probabilistic interpretation of the DIS cross section is thus precluded. 

Our paper thus explains the origins of nuclear shadowing and leading-twist diffrac- 
tion, giving a new, first principle, perspective on these problems. Our formalism 
of final-state interactions has recently been used to analyze single-spin asym- 
metries in deep inelastic processes and to show that such asymmetries survive 
in the Bjorken limit, contrary to conventional arguments which claim that final 
state interactions are always power-law suppressed in the large scale of hard QCD 
processes 0. 



2. The Foundations of the QCD-Improved Parton Model 

Soon after the observation of Bjorken scaling (and before the advent of QCD) 
it was suggested 0] that the DIS cross section is fully determined by the target 
wave function. Specifically, consider the Fock expansion of the nucleon state 
|A^) in terms of its quark and gluon constituents at equal Light-Cone (LC) time 
T = t + z/ c = = y'^ + , 



\N) 



n 



i'uudixi, k_i_i, Xi)\uud) 



+ ^wg(- • ■)\uudg) + ... + . .)\uu. 



+ 



(2) 



Each Fock state \uud . . .) is weighted by an amplitude if) which depends on the LC 
momentum fractions Xj = kf /p^ ijli^i = 1)? the relative transverse momenta 
ku k±i = 0), and the helicities Aj of its const it uentsQ. The DIS cross section 
thus appeared to measure the single parton probabilities Vj/N{xB, Q^) as defined 
in (|l|), which express the probability for finding (at resolution 1/Q) a parton j 
carrying the momentum fraction Xb = Q'^/'^P ■ Q of the nucleon. Here q is the 
virtual photon momentum (g^ = — Q^) and p the target nucleon momentum. 

Later analyses of perturbative QCD (PQCD) have established the QCD fac- 
torization theorem to all orders in the coupling. The DIS cross section can be 
expressed for each parton type as a convolution of a perturbatively calculable hard 
subprocess cross section and a target parton distribution. The parton distribu- 
tions are given by operator matrix elements of the target. For the (spin- averaged) 

"'^See Ref. for the normalization conventions. 
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quark distribution in the nucleon of momentum p, 



fq/N{xB,Q^) = ^ dy exp{-ixBP^y ) 




(3) 



I 

Jo 



X 



{N{p)\q{y )7+Pexp ig 



dw~A^{w') q(0)|A^(p)) 



where all fields are evaluated at equal LC time y"*" = and vanishing transverse 
separation = 0. The light-like distance between the absorption and emission 
vertices of the virtual photon in the forward amplitude is measured by y~ . The 
path-ordering P orders the gauge fields according to their position on the light- 
cone and ensures the gauge invariance of the matrix element. 

The identification of the quark distribution (H) as a probability distribution (|l|) 
is made in LC gauge n ■ A = = 0, where the path-ordered exponential in (|^) 
reduces to unity, and one finds /q/iv V^/n- A recent derivation in the more 
general case of non-forward matrix elements (Skewed Parton Distributions) may 
be found in Ref. [0. Thus the DIS cross section appears to be fully determined 
by the probability distribution of partons in the target. However, as we shall 
show the expression for f^^/N cannot be given by (|) in LC gauge. 

In a general gauge the matrix element (|^) depends on final state interactions (FSI) 
of the struck quark with the gauge field of the target via the A^-dependence of 
the path-ordered exponential. Based on the above argument in LC gauge, it 
is generally believed that the exponential is a gauge artifact and thus that the 
presence of FSI does not infiuence the cross section. But this assumes that fq/N 
is given by (|^) in all gauges, including LC gauge. Here we find that final state 
rescattering in fact does change the DIS cross section in all gauges. Our analysis 
is consistent with the QCD factorization theorem and with the form of the 
parton distributions in all gauges except LC gauge. 

The infiuence of FSI we find at leading-twist is specific to gauge theories. The 
impossibility to interpret parton distributions as probabilities could thus not 
be inferred before the advent of QCD. Instead, the equivalence between DIS 
structure functions and the target wave function was assumed, though it was 
only shown in a theory with Yukawa coupling [|l|. 

The expression for /q/^v is valid for covariant gauges in the Bjorken limit, 
which selects the A'^ field of the target. We shall show that setting then A'^ = 
in d^) leads to an incorrect expression for fq/N- From a mathematical point of 
view this means that the high energy Bjorken limit does not commute with the 
A"*" — > limit of LC gauge. In fact (see section 7) the high energy and LC gauge 
limits do not commute even for ordinary elastic electron scattering. 

In section 3 we recall why in Feynman gauge final state interactions among the 
spectator partons of the target system do not affect the DIS cross section at lead- 
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ing twist. We then show that this general argument does not apply to rescattering 
of the struck quark. 



In section 4 we discuss the Glauber-Gribov picture and show why it implies that 
the final state interactions, resummed in covariant gauges by the path ordered 
exponential of @, affect the cross section. We then study a simple perturbative 
model of rescattering effects in section 5, for which explicit expressions of the 
amplitudes can be obtained at small xb- Using this example we demonstrate in 
section 6 that rescattering of the struck quark on the target can cause a leading 
twist shadowing effect. 

The analysis of sections 3 to 6 is carried out in Feynman gauge. In section 7 we 
show why rescattering effects can persist even in = gauge, in contradiction 
with the form of the matrix element. As is well-known, this gauge is singular 
- in particular, the gluon propagator 



k"^ + ie 



^ n ■ k 



(4) 



has a pole at n-k = k~^ = which requires an analytic prescription. In final-state 
elastic scattering of the struck quark the exchanged momentum k~^ is of O {1 / u) 
in the target rest frame, which enhances the second term in the propagator (^. 
This enhancement allows rescattering to contribute at leading twist in LC gauge. 

We reevaluate our model amplitudes using LC gauge in the Appendix. Although 
the expressions for the individual diagrams depend on the prescription used at 
n ■ k = 0, the prescription dependence vanishes when all diagrams are added. 
The scattering amplitudes which we calculate up to two-loops in LC gauge thus 
agree with the result in Feynman gauge. 

For the issues of this paper, the spin and color of the quarks are not relevant. 
We therefore conduct our discussion in the simpler framework of abelian gauge 
theory with scalar quarks. 



3. Effects of final state interactions in deep inelastic scattering 

The DIS cross section is given by the discontinuity of the forward amplitude, 

a(7*r ^ X) = ^Disc^ (7*^ ^ ^T) (5) 

where M is the target mass and u the photon energy in the rest system of the 
target. We take the Bjorken limit I'yQ^ = —q'^ oo with xb = Q'^/2Mv fixed. 
In the LC notation k = (/c+, k~, k±), where k^ = k^ ± k^, the photon and target 
momenta are (at leading order) 

q = {-MxB, 2u,0±) 
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p = (M,M,Ox) 



(6) 



In the following we define a final state interaction (FSI) as any interaction which 
occurs after the virtual photon has been absorbed. Here 'after' refers to LC 
time, = + y^, in the frame (^. In deep inelastic scattering initial state 
interactions (ISI) occur only within the target bound state and determine the 
target wave function (|^). We shall show that soft rescattering of the struck 
quark in the target also affects the DIS cross section. 

We can distinguish FSI from ISI using LC time-ordered perturbation theory 
(LCPT) Fig. 1 illustrates two LCPT diagrams which contribute to the 

forward 7*T — > 7*T amplitude, where the target T is taken to be a single quark. 
We use these diagrams in a generic sense here, while in sections 5 and 6 we con- 
sider them in the framework of a specific perturbative model of the DIS process. 




Figure 1: Two types of final state interactions, (a) Scattering of the antiquark 
{p2 line), which in the aligned jet kinematics is part of the target dynamics, (b) 
Scattering of the current quark {pi line). For each LC time-ordered diagram, 
the potentially on-shell intermediate states corresponding to the denominators 
Da, Dfj, Dc are denoted by dashed hues. 

We recall that in LCPT the ' — ' momentum component is not an independent 
variable, but is given by the on-shell condition, k~ = {k\ + m?)/k^. Each 
propagating line has a factor l/Zc'*', and there is a denominator factor 

Dint = Y.k- -Y.^- +ie (7) 

inc int 

for each intermediate state, which measures the LC energy difference between 
the incoming and intermediate states. In Feynman gauge (which we use in this 
section) an imaginary part or discontinuity can arise only via the ie prescription 
in (^, when LC energy is conserved and the intermediate state is on-shell. 
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We consider the 'aligned jet' (or parton model) configuration |T^, where the hard 
vertex is taken at zeroth order in the strong coupling: 7*q q. In the aligned 
jet kinematics the momentum pi of the struck quark in Fig. 1 is the only one 
which grows in the Bjorken limit: pj~ ~ 2z/, with pi± independent of u. All 
momenta in Fig. 1 other than q and pi remain finite in the Bjorken limit. The 
condition that the momentum fraction of the struck quark equals Xb follows from 
the conservation of '+' momentum, given that pf = 0{l/v). 



We recall (see, e.g., Eq. (A5) of Ref. ||T3|) that the virtual photon polarization 

vectors may be chosen as 

e{\ = ±l) = --L(0,0,l,±z) 

£(A = 0) = ^(1,-1,0,0) (8) 

Since we take all lines (except the gauge bosons) in Fig. 1 to be scalars, the 
longitudinal photon coupling e{X = 0) ■ {pi + ki — P2) — Q dominates over the 
transverse ones in the Bjorken limit. The two longitudinal photon couplings 
together contribute a factor to the forward amplitudes in Figs, la and lb. 

Both diagrams in Fig. 1 contain final state interactions between the 7* vertices. 
Only the three intermediate states indicated by dashed vertical lines can kinemat- 
ically be on-shell and thus contribute to the discontinuity of the diagrams via the 
vanishing of the corresponding denominator Da, or D^. We wish to ascertain 
whether the sum of these discontinuities gives a leading-twist contribution to the 
DIS cross section through the optical theorem (^. We use Feynman gauge in the 
following discussion. As we shall see in section 7 and Appendix C, the specific 
Feynman diagrams causing FSI effects in DIS actually depend on the gauge. 

The three denominators of Fig. la are 

Da = -Pi -P2 -{p-h)~ 

= 2u- + +M- + - + 
^ pt pt M-kf 

(9) 

pI^ + Pi± + kl^ (fci± + k2A.f + 

Uh — /z/ x \- ivi x m 7T rx m — 

pt pt kt M-kt-kt 

_ pI± + {p2i_ + fc2±)' + _ (fci± + k2±? + 

pt pt + kt M-kt-kt 

and have the form 

Da,b,c = 2U - + fa,b,c (10) 

Pi 
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where fa, fb, fc are independent of u in the ahgned jet configuration. If we consider 
these denominators as functions of pf then the three conditions -Da,fe,c = give 
to leading order the same value of pf, 



Pi 



pIj_ + 171^ 

2^ 



Oil 



(11) 



All denominators and other factors in the LCPT expression of Fig. la except 
Da, Df, and Dc are insensitive (at leading order) to a relative change in pf of 
O [l/u). Thus, as far as the discontinuity of Fig. la is concerned, we can regard 
the other factors as constants in the p^^-integral containing the denominator poles, 

Disc(Fig. la) oc Q^Disc / ^^77^ . . (12) 

where the factor stems from the photon couplings. All remaining factors in 
the proportionality are independent of u. Each of the three denominators in (0) 
gives a //-independent contribution to the discontinuity in the Bjorken limit. This 
means that each partial discontinuity contributes to the DIS cross section of Eq. 
(^ at the leading twist level, cr{j*T ^ X) oc l/Q^. However, as is easily seen, 
the contributions from the three poles which to leading order occur at the same 



value (11) of pi cancel at leading twist. 



The above argument is generic and applies to arbitrarily complex diagrams hav- 
ing no interactions on the current quark line pi. The remarkable fact that FSI 
between target spectators do not affect the DIS cross section only relies on the 
Bjorken limit, which as u 00 provides an 'infinite energy reservoir' which 
compensates any target excitations. 

The situation is quite different for diagrams like Fig. lb where the current quark 
reinteracts. In (quasi-)elastic scattering of the current quark the momentum 
transfer fc^ oc 1/z/. We may check explicitly that this range of momentum trans- 
fer indeed gives a leading-twist contribution to each partial discontinuity. The 
denominators are now of the form 

P{ 

D, 2.-?i±-t^-|i + 9. (13) 

Pi 1^2 

" 2. + g. 

where ga,h,c are again independent of z/. For example, picking up the Da = Q pole 
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in the pf integral we have 

where pi is given by ([llD and the factor pi originates from the interaction in 
Feynman gauge. 



Note that and Df. are still of O {v) at the value (|TT|) of Pi for which Da = 0. 
The fact that the contributions from Da = 0, Db = and Dc = thus occur at 
distinct values of pf means that they no longer cancel. DisCa is independent of u 
and hence contributes to the DIS cross section at leading twist. We conclude that 
rescattering of the current quark generally affects the cross section. In section 
6 we shall demonstrate, in terms of an explicit perturbative example, that this 
conclusion is indeed correct. 

Since the LC energy differences D{,^c ocu at Da = 0, the struck quark rescattering 
occurs on the light-cone, ~ O (l/u). This rescattering is part of the dynamics 
described by the path-ordered exponential in the matrix element @, where all 
A'^ fields are evaluated at the same LC time y~^. During its passage through the 
target the struck quark has no time to emit or absorb gluons, it only 'samples' the 
Coulomb field of the target. The rescattering nevertheless changes the transverse 
momentum of the quark and influences the cross section. This is analogous to the 
LPM effect which suppresses the bremsstrahlung of a high energy electron 



in matter due to Coulomb rescattering within the formation time of the radiated 
photons. 



4. The Glauber-Gribov Picture of Shadowing 

DIS data on nuclear targets A has shown that nuclear structure functions are 
suppressed for xb ^ 0.05: F^{xb,Q'^) < AF^{xb,Q'^) This is generally 

interpreted as a leading twist 'shadowing' effect, arising from quantum mechanical 



interference 0, |T0| . The coherence length of the virtual photon in the target rest 

{y-) (15) 



frame (g) is long at small xb, 

2z/ 



g2 Mxb 

as can also be seen from Eq. (^. Rescattering from different nucleons in the 
nucleus can thus interfere. 

In the aligned jet kinematics the virtual photon fluctuates into a qq pair with 
limited transverse momentum, and the (struck) quark takes nearly all the longi- 
tudinal momentum of the photon. Using the notation of Fig. 1, where the initial 
q and q momenta are denoted pi and P2 — ki, respectively, we have 

p^ ~ 2u 
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pt -K ~ -Mxb (16) 
Pi± = -{P2j_ -hj ~ Aqcd 

The (covariant) virtualities pi and (p2 ~ ^i)^ are limited. Hence {pi +P2 — ~ 
Piipt ~ kf) ~ — as required by momentum conservation. The virtual quark 
pair is put on-shell by a (total) momentum transfer k from the target, with 

k+ = pt + pt - ^ pt + Mxb (17) 

The DIS cross section is dominated by minimal transfers k~^, which for the final 
antiquark momentum gives 

pt ~ Mxb (18) 

With this kinematics in mind the Glauber-Gribov picture of shadowing can be 
summarized as follows. At small xb the antiquark momentum p2 oc 1/xb is large 
but the momentum transfer ~ Mxb is small. The scattering will therefore 
have a diffractive component. In particular, the quark pair may scatter elastically 
on a 'front' nucleon A^^i in the nucleus before suffering an inelastic collision at a 
'back' nucleon N2, as indicated on the Ihs of Fig. 2. The small momentum 
transfer k~^ at A^^i required to put the quark pair on-shell can be absorbed by 
the nuclear wave function. Hence this amplitude interferes with the amplitude 
for a single scattering on N2 shown on the rhs of Fig. 2. The interference is 
destructive due to the imaginary nature of the Pomeron exchange amplitude at 
Ni and the factor of i resulting from the intermediate state between Ni and N2 
going on-shell. 




Figure 2: Glauber-Gribov shadowing involves interference between rescattering 
amplitudes. 

This shadowing effect on the DIS cross section is not compatible with the cross 
section being determined by the parton probabilities V of Eq. (0). Since the 
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Pomeron amplitude in Fig. 2 is imaginary it must involve on-shell intermediate 
states. But initial state interactions in the target before the virtual photon is 
absorbed cannot create on-shell intermediate states - they would constitute decay 
channels of the target. We conclude that Glauber-Gribov shadowing involves final 
state interactions and hence must be associated with the path ordered exponential 
in (i. 

5. A Perturbative Example of Shadowing 

We shall construct a perturbative example of the physics of Glauber-Gribov shad- 
owing, which is simple enough to allow explicit expressions for the scattering 
amplitudes at small xb- We use this example in section 6 to verify the general 
result of section 3 that final state interactions between target spectators do not 
affect the DIS cross section, whereas rescattering of the struck quark does. 

In this section we use standard covariant perturbation theory in Feynman gauge 
of a scalar abelian gauge theory. 

We consider the forward 7*T 7*T amplitude of Fig. 3, the discontinuity of 
which gives a contribution to cr(7*T) at order aaf via the optical theorem (^. 
Since we may assume the charges of the target T and the 'quark' q to be distinct, 
we can focus on the gauge invariant set of diagrams in which the gluons are 
exchanged between the quark pair and the target. Each gluon ki can couple to 
either the q or the q line, and all distinct permutations of the gluon vertices are 
included. 




T(p) T(p) 

Figure 3: Forward 7*T — > 7*T amplitude. All attachments of the exchanged 
gluons to the upper scalar loop are included, as well as topologically distinct 
permutations of the lower vertices on the target line. 

Taking the discontinuity between gluons and k^ gives a contribution which 
models the interference term of Fig. 2. The scattering on N2 is given by single 
gluon exchange, while the Pomeron exchange on Ni is modelled by two gluon 
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exchange. The discontinuity between gluons k2 and gives the square of the 
'Pomeron' exchange amphtude. We calculate the one-, two- and three-gluon 
exchange amplitudes for 7*T — > qqT explicitly for <C 1, making use of the 
results of Ref. fl^ where a similar model was studied. Since the target T is taken 
to be elementary this model does not have shadowing in the conventional sense 
described in section 4. It nevertheless demonstrates how final state interference 
effects reduce the DIS cross section. 

We work in the target rest frame @ and in the aligned jet kinematics of Eqs. 
(|16D and (|18]). The Feynman gauge calculation is simplified by assuming]^ a large 
target mass M. Hence the kinematic limit we consider is 

2z/ ~ p7 > M > > ki±, p2±, k', m > kf , /c+ = Mxb + pt (19) 

where m is the mass of the q, q quarks and k = J2iki is the total momentum 
transfer from the target. 



5.1 Single Gluon Exchange Amplitude A 



The three Feynman diagrams are shown in Fig. 4. As in section 3 we use the 
virtual photon polarization vectors (|D and find that the dominant (leading twist) 
contribution comes from e{\ = 0)-pi — Q. Diagram 4c is proportional to e-{p+p') 
and is thus non-leading. Diagram 4a involves the quark propagator 

(P2 - kf - ~ p-{p+ - k+) - (p2± - k^f -m^ = -D{p2± - k±) (20) 

where we used (0) and defined 

D{p^) =P^Mxb+P±'^ + m'^ (21) 

Similarly the quark propagator in diagram 4b gives D{p2±). The full amplitude 
in the limit (|T9t) is 



MP2 ,P2±,k±) 



2eg^MQp^ 

k^. 



1 



1 



^(m) D{p2i. - ki) 



(22) 



We may readily verify that this contribution is of leading twist. The l+T f!+X 
DIS cross section is[|, 

^4 da _ a 1-y 1 r dp^ d'^p2± d'^k± ^ 

^ dQ^dxB 167r2 2MuJ p^ (27r)2 (27r)2 ' ' ^ ' 

^The expressions for the scattering amplitudes that we derive at large M are actually valid 
also when M and k± are of the same order. This is seen directly for the Born amplitude of 
Fig. 4, and from the LC gauge calculations in the Appendix for the loop amplitudes. 

•^Here the lepton £ is assumed to have spin ^ . 
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y*(q) 





(a) 



(b) 



(c) 



Figure 4: Single gluon exchange diagrams in scalar abelian theory. 

where y = v / E^. The factor in \A\'^ combines with l/2Mv in (^) to make 
the rhs independent of in the Bjorken limit, when the soft momenta k± and 
P2 are integrated over any finite domain. 

We also note that the dominant contribution to the DIS cross section at small xb 
comes from ~ Mxb and P2 ~ (P2± + im?)/MxB as assumed in (|18|). To see 
this, note that the amplitude A oc P2 for p2 ^ {p2j_ + m'^)/MxB, while A oc l/p2 
for P2 > {pIj_ + m?)/MxB- 

Since A oc l/fcj, for /c^ — the cross section ( PHj ) has a logarithmic singularity in 
this limit, which is regulated by the longitudinal momentum exchange at kx^ ~ 
fc+ ~ Mxb- This logarithmic behavior occurs only at lowest order and will 
not be relevant for our conclusions. 

It is instructive to express the cross section also as an integral over the transverse 
distances r±, R± conjugate to p2±, k±. Defining 



A{p^,r^,Rj 
and using 



cPp2± (fkj 



(27r)2 (27r): 



^(P2 5 P2±, k±) exp (if± ■ p2± + iR± ■ (24) 



V{mr±] 



(27r)2p^ + m2 27r 
where Kq is a Bessel function, and 



(2vr) 



k'i 



1 

2^ 



log 



\R± + r±\ 
Ri 



we get from (|2 



where 



mn = P2MXB + 



(25) 

(26) 

(27) 
(28) 
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The contribution (^) to the DIS cross section can then be expressed as 



da 
dQ"^ dxB 



a^al I — y 



xbM"^ j dp2 



P2 



(j)2 Mx 



B + m' 



X 



d'^u±d'^U^ 



Ko{u±)\og 



\U± + u±\ 



(29) 



Here the dimensionless integration variables were defined as u± = f±m\\ and 
U± = R±m\\, showing that the typical transverse distances f±, R± scale as l/m\\. 
The P2 integral in (|29|) is logarithmic^ at large > m? /Mxb-, where the aligned 
jet 7*q q subprocess turns into 7*7 — > qq [|T3|. 



5.2 Two-Gluon Exchange Amplitude B Fig. 5 shows two of the altogether six 
two-gluon exchange diagrams which give leading contributions to the 7*T — qqT 
amplitude for -C 1 in Feynman gauge. Diagrams with 4-point vertices (c/. 
Fig. 4c) are again suppressed in this gauge. We illustrate the calculation of this 
one- loop amplitude using the diagrams of Fig. 5. 




Figure 5: Double gluon exchange diagrams. In Feynman gauge four more dia- 
grams contribute at leading order, where one or both of the exchanged gluons 
attach to the quark (pi) line. 

Our assumption (^) of a large target mass M simplifies the loop integral by 
suppressing the momentum components. For the overall exchange we find 
from the mass-shell condition p'"^ = {p — kY = that 

k' = kl + kl = -^<^k^, k^ (30) 

The corresponding suppression for the loop momentum k^ ~ — fcg results from 
the sum of the uncrossed and crossed gluon attachments to the target line in Fig. 

^ We also note that ( p9| ) contains a coUinear singularity when to ^ 0. In this limit the 
exchanged gluon becomes a collinear line in the language of Ref . Isj . 



13 



5, 



-ig2Mf 



{p - kif -M'^ + ie^ {p- -M^ + ie 



2ig^M 



1 1 

+ 



ki — ie k^ — ie 



2ig^M 2m5{k^^] 



Making use of Eqs. (pof ) and (|3T|) we find 

2eg^MQp2 r d%±dkt 1 



(31) 



B{5a) + B{5b) 



Dip 



'2± 



ki 



(27r) 



1^1 1^2 



1 



1 (32) 

H - (2p2± ■ k2± - kl)/p2 - ie 

Symmetrizing the integrand in ki ^ k2 and recalling (|l^) the last factor becomes 



Thus 



H - (2P2± ■ k2± - kl)/p2 - ie k^ - k+ + {2p2± ■ - kj)/p2 + ie 
~ i7r(5(A;J) 

ieg'^MQ P2 



Biha) + 5(56) 



2± 



D{p2i. - k^) J (27r)2 kl^ k. 



(33) 



(34) 



2± 



Adding the contributions from the remaining four diagrams we find for the full 
two-gluon exchange amplitude 



B{P2 ,P2±,kj 



Wmqp2 
1 



cPk 



1± 



(2vr)2 kl^kl^ 
1 



(35) 



_D{p2i_) D{p2i_ - k 



1±J 



D[p 



'2± 



k 



+ 



2±j 



2± 



where fc2± = k^ — k^^. We note that the amplitude is fully imaginary as required 
by crossing symmetry, since i? oc P2 — >■ cxo and two gluon exchange has even 
charge conjugation. Thus our model captures the essential features of Pomeron 
exchange. We note also that B{p2 ,P2±, k±) oc logfcj, for k± 0. In contrast 
to the single gluon exchange contribution to the DIS cross section, the square of 
(|35[) can thus be safely integrated over k± and (for m ^ 0) over • 



Due to conservation of the transverse distances f±, in the peripheral scatter- 
ing, the Fourier transform (E^) returns the simple form 

^2 



B{p2,r±,Rj 



-ieg*MQp2V{m\\r^)W\f^, R 



2! 



WA 



(36) 
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where we used (^) and (|26|). 



We stress that in the — > hmit, the amphtude B is dominated by the config- 
uration where the intermediate state between the two exchanges is on-shell. This 
can be seen by calculating B in LC time-ordered perturbation theory, where this 
intermediate state is associated with a vanishing denominator (^. Alternatively, 
we may note that since the real part of B is suppressed in the xb —>■ limit the 
full amplitude is (via the optical theorem) given by its discontinuity. This is true 
in all gauges since B is gauge invariant. 



5.3 Three- Gluon Exchange Amplitude C 

No qualitatively new aspects appear in the calculation of this two-loop amplitude. 
Permuting the attachments of the three gluons on the target line one finds in 
analogy to (|3l|) that ~ for all exchanges {i = 1,2,3). Similarly the kf 
integrations are simply evaluated after symmetrizations analogous to (p3D. The 
final expression in momentum space is 

C{p. ,m, k,) = --eg MQp, J 

1 3 3 
X — \ — 

D{p2±) D{p2^ - h^) D{p2^ - fcix - k2±) 

where /c3j_ = k_i — ki^ — k2±. 
The Fourier transform ( ^4|) gives the amplitude in transverse coordinate space as 

Cip2, rl, R±) = -yg''MQp^Vi7rn\r^)W%f^, R^) = t^w^A (38) 



(37) 

1 

D(p2±-k^). 



Similarly to the B amplitude, C arises from the intermediate states between the 
rescatterings being on-shell in the — > limit. Again this must hold also in LC 
gauge. Since in the two-loop case there are two consequtive intermediate states, 
C is purely real. 

), it is apparent that the sum of gluon- 



From the expressions (|27|) , ( ^6|) and 
exchange amplitudes exponentiates. 



M{p2,f±,R±)=A + B + C... = -2ieMQp:^V 1 - exp{-ig^W) (39) 



As noted at the beginning of this section, we have assumed the charges of the 
quark and target lines to be distinct. This allows us to restrict our analysis to 
the subclass of Feynman diagrams considered above, since diagrams with different 
powers of the charges cannot cancel in the DIS cross section. However, we should 
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note that at the level of three gluon exchanges there are new types of diagrams 
which have the same charge dependence as C in Eq. (p7|). For example, one of 
the three gluons may be exchanged between the quarks while another forms a 
loop on the target line. The fc^-dependence of this contribution would differ from 
that of (^). We do not further consider such contributions. 



6. Effects of Rescattering on the DIS Cross Section 



We now use our perturbative amplitudes to demonstrate that final-state rescat- 
tering of the struck quark affects the DIS cross section. In the previous section 
we used covariant (rather than time-ordered) perturbation theory, and thus did 
not distinguish between initial (ISI) and final (FSI) state interactions. However, 
diagrams involving rescattering of the struck quark necessarily are FSI because 
the exchanged gluon couples to the struck quark {pi) line after the virtual photon. 
We shall see that precisely such diagrams contribute to the cross section. 



We consider the DIS cross section 
distances r±,R± defined in (El), 



3D expressed as a sum over the transverse 



da 



a I — y 1 
dQ^ dxB ~ 167r2 2Mv 



^d'f^d'R^\M\' 

P2 



where 



\MiP2,r±,RA 



sm 



is the resummed amplitude 
respectively. 



MP2,'^±,R2 

g'Wif^,R^)/2 

and V, W are given in Eqs. 



(40) 



(41) 



51) and 



The fact that the coefficient of A in (^Tj) is less than unity for all r±, R± shows 
that the rescattering corrections included in Ai reduce the cross section. This 
effect agrees with the Glauber-Gribov picture of DIS shadowing and must be 
present also in LC gauge (see section 7). 

The forward 7*T '~f*T amplitude in Fig. 3 can also be cut through some of 
the gluon lines, corresponding to final states with real gluons. Such contributions 
have, however, a different target mass M dependence (c/. Eq. (|30|)). Similar 
arguments suggest that other contributions, even if they are of the same order in 
the coupling constants, cannot change the conclusion that the DIS cross section 
is influenced by final state interactions. 

In section 3 we gave a general argument (in Feynman gauge) which showed that 
final state interactions between target spectators cannot influence the DIS cross 
section (c/. Fig. la). We shall now check this statement using our perturbative 
amplitudes. 
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In the aligned jet kinematics the antiquark belongs to the target system. We thus 
consider the subset of diagrams like Figs. 4a and 5 where all exchanged gluons 
attach to the q (^2) line. One can easily verify that this subset of diagrams is 
gauge invariant in the class of covariant gauges in our kinematic limit (plQl) . The 
corresponding sum of cuts in Fig. 3 is then proportional to 



|5q|2 + 2Re(Aqq) 



(42) 



where the subscript q indicates the subset of diagrams. 



Diagrams where all gluons attach to the antiquark line can involve both ISI 
and FSI. Since the two-gluon exchange contribution ([3^ ) is imaginary it must, 
however, involve rescattering of on-shell intermediate states which can only arise 
after the virtual photon has been absorbed. Similarly the (real) three-gluon 
exchange amplitude C (|37|) involves double rescattering of on-shell states. Hence 
all our amplitudes (except the Born term A) involve FSI. 



q, i?q, Cq coutributious to the expressions 



It is straightforward to identify the A, 

(p2D, (0), (|37|) of the full one-, two- and three-gluon exchange amplitudes in 
momentum space. According to Eq. (^) the antiquark propagator next to 
the virtual photon vertex gives a denominator D{p2± — k±) for all diagrams in 
our subset. This factor appears explicitly in each amplitude. Dimensionally 
regularizing the logarithmic infrared divergencies at ki±_ = we thus find 



Sq{P2 ,P2±,k 

where 

R2{k±) = 



eg^MQp-2 



{2TT)W{p2±-k^)\ 



12 4 



[R2{k^)Y - 3^i3(A;j 



(43) 



TC 



D/2 



r(|-i)jV(2-f) 



4-_D 



T{D-2) 



(44) 



Ri3{kj 



1 



d^ki^ d^k2± 






^T{3-D) 


khklAk±-h±-k2±)^ /cl"'^ 




r(f 


-3) 



Expanding R2 and -R13 around D = 2 gives 

r 4 



R2ik, 



7C 



{ + + 1^ (^^' " - 2) + 



+ 24 



A;j-^ [D-2 
1 



2/ _^7r^ _ 14^(2)(i) {D-2Y + O {D-2) 



(45) 



17 







-2D 


1 


4^3 



12 



iD-2) 



127 
D-2 



TT 



+ 67 - — + 



4^3 _ 16^(2) (1) (D-2) + C (£'-2)' 



where 7 ^ 0.577 is Euler's constant and ip^^\z) is the {n + l)*'^ logarithmic 
derivative of the gamma function with ■ip^'^^l) ^ —2.40. Hence 



o k 



D 



8-2D 



6^(2) (1) {D-2) + 0\{D-2f 



(46) 



vanishes at D = 2. Thus FSI between the target spectators do not change the 
DIS cross section. 

We conclude that in covariant gauges, only final state interactions which involve 
rescatterings of the current quark affect the DIS cross section. 



7. Light-Cone Gauge = 

We have seen that in a covariant gauge, the DIS cross section is influenced by 
final-state interactions of the struck quark in the target. This soft physics is 
contained in the path-ordered exponential of the matrix element (^) in a general 
gauge and appears to vanish in LC gauge, n ■ A = = 0. 

However, as we have seen in section 5 the amplitudes B and C arise from on- 



shell intermediate states in the limit. Thus in (^Of ) the contribution 

of |A^P, whose expansion starts as |-Bp + 2AC, also arises purely from on-shell 
intermediate states. The presence of such on-shell states is gauge independent 
and they can only occur in the final state. We conclude that the DIS cross section 
is influenced by final state interactions in all gauges. Thus parton distributions 
cannot be fully determined by parton probabilities in the target. 

Let us now discuss some features inherent to the LC gauge preventing parton 
distributions from being probabilities, in other words making the expression (|^) 
for /q/AT incorrect in LC gauge. It turns out that terms which are next-to-leading 
corrections in a general gauge cannot be ignored in LC gauge. To see this, it is 
helpful to recall how the exponential arises from perturbative diagrams. 

As explained in Ref. each quark field is associated with an ordered exponential 



[A+] = P exp 



00 



ig / dy A+{y' 







(47) 



where the gauge field ^4+ is evaluated on the light-cone, y~^ = y± = 0. This factor 
arises from the interactions of the struck quark as it moves through the target. 
While the path in ( ^Tj ) extends to infinity, there is a partial cancellation between 
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the quark fields in the matrix element (0) leaving a path of length ~ 1/Mxb 
equal to the coherence length of the virtual photon. Only interactions within this 
LC distance can influence the cross section. 

Expanding the exponential gives 
[A+] = l + ig dy^A+ivY) I + ig dy^A+{y^) + .. 

Jo I Jy^ 

(48) 

«^ dkt A-^{kt) 2 r dkfdkt A+{kt)A+{kp 



+ 



where 



dk{ ) 2 f 

J-oo 27r kt -le^^ J-oo (27r)2 (kf + k^ - ie){kt - te) 

A+ik-')expi-tk+y-) (49) 

-oo ZTT 



Figure 6: Scattering of the struck quark on the gauge field of the target which 
gives rise to the ordered exponential (^). 



The terms in the expansion ( ^q) arise from the perturbative diagrams of Fig. 6, 
where the cross indicates the virtual photon vertex. The struck quark momentum 
is asymptotically large, Pi —>■ oo, implying that the quark moves along the light- 
cone y~^ = y± = 0. The two-gluon exchange term in Fig. 6 is given by 

,...^2..2 piA^kt)piA^kt) 2 A^kt)A^kt) 



Vi{kt + k+) + ie] i-pik+ + le) + kt - ie){kt - is) 

Thus we find equivalence to the expression (^) by approximating {2pi — ^2) • 
A{k2) — Pi A'^ {k2) 1 i.e., by keeping only the asymptotically large component of 
Pi. This is correct in all gauges except A^ = 0, where this 'leading' term actually 
vanishes. 

Neglecting the dependence of the matrix element (^ on the gauge field A{k~^) 
in LC gauge is equivalent to assuming that interactions of the struck quark with 

the gauge field such as {2pi — k2)± ■ A± do not contribute at leading twist. The 
following example shows how this assumption can fail. 
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As a simple illustration of how the high energy and LC gauge limit can fail to 
commute we consider the elastic process q(pi — k) T{p) q(pi) T{p — k), where 
p = (M, M, 0_l) and ^ at fixed pi±, k±. Momentum conservation implies 

^+ ^ {2pi± -k^)-k± ^ Q 
Pi 

k- = --^ fixed (51) 

The interaction of the gauge field with the quark is given by {—ig){2pi — k)^ ■ 
d^^{k). In Feynman gauge the propagator is d^{k) = —ig^^/{k'^ + ie) and the 
coupling is dominated by —igpidp^{k), which is analogous to the interaction 
in the ordered exponential. The elastic amplitude 

A{qT ^ qT) = -2g^M^ (52) 

is thus oc Pi as befits Coulomb exchange. 

In LC gauge the propagator (^ satisfies d~j^Q{k) = 0, hence the p~[ component 
does not contribute. Yet the elastic amplitude is gauge independent and must 
still be given by (|52|) . The absence of the factor in the numerator coupling 
is in fact compensated by the factor k~^ oc 1/pi in the denominator of the LC 
gauge propagator (|p. The dominant contribution is from — (2pi — k)^ ■ d~^{k) 
and the result indeed agrees with (|5^). 



Note that if we had kept only the d~^'^{k) part of the gauge propagator in the high 
energy limit and then chosen LC gauge the elastic scattering amplitude would 
have seemed to vanish. This incorrect result is analogous to the apparent absence 
of rescattering effects in the matrix element (^ for A'^ = 0. 

In the Feynman gauge calculation of section 5 we saw that the reinteractions of 
the struck quark with the target are essentially elastic, the intermediate states 
being on-shell in the xb ^ limit. It is thus not surprising that the calculation 
of the scattering amplitudes in LC gauge has many features in common with 
the elastic scattering example above. Details of the calculation of the one-loop 



and two-loop amplitudes B and C (^) and (^) in LC gauge are given in the 
Appendices. 

In LC gauge the Feynman rules must be supplemented with a prescription for the 
A;^ = pole of the propagator (1). Three prescriptions that have been studied in 



the literature |TT], |1^, are given in Eq. (^) of Appendix A. The contributions 
of the individual diagrams shown in Fig. 7 for the one-loop amplitude B depend 
on the prescription. However, the k^' = poles cancel when all diagrams are 
added. Their sum is thus prescription independent and agrees with the Feyn- 



man gauge result (PH). We verify the prescription independence of the two- loop 
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amplitude C in Appendix B. A consistent procedure for regulating the spurious 
poles is also discussed there. 

As we have already emphasized, final state interactions (FSI) modify the DIS 
cross section also in LC gauge due to the presence of on-shell intermediate states 
between the rescatterings in the amplitudes B and C . However, while in Feyn- 
man gauge it is the rescattering of the struck quark pi which affects the cross 
section, in LC gauge those rescatterings actually do not contribute. Indeed, we 
show in Appendix C that contributions from diagrams like Fig. 7c and 8b to the 
individual amplitudes cancel in the cross section. Thus in LC gauge, indepen- 
dently of the prescription, the cross section is modified by FSI occurring on the 
antiquark p2i ^-c, within the target system. Choosing the = gauge shifts 
the rescatterings of the quark present in Feynman gauge to rescatterings of the 
antiquark. As also shown in Appendix C, in LC gauge the partial amplitude 
where only attachments to p2 are kept equals the full amplitude, up to a phase 
factor. Which particle actually scatters in the quark-antiquark system depends 
on the gauge, but the presence of on-shell intermediate states does not. 

Subtleties can appear when using the Kovchegov (K) prescription (see Eq. (|65D), 
since the imaginary part arising from a physical cut can be changed by the imag- 
inary part created by the prescription itself. The K prescription simulates the 
physics of the rescattering corrections by introducing an external gauge field into 
the dynamics. Unlike the Principal Value (PV) or Mandelstam-Leibbrandt pre- 
scription, the K prescription is not causal, and thus it would normally not be 
used for solving the bound state problem and light-cone wave functions of an 
isolated hadron in QCD. The solutions for the light-cone wave function of the 
target hadron in the presence of an external gauge field can have complex phases. 
This is apparently the way in which the light-cone wave functions of a nucleus in 
the Kovchegov light-cone gauge prescription mimic the effects of rescattering of 
the fast quark and the Glauber-Gribov shadowing modifications of the structure 
functions. If this picture could be validated, the Kovchegov LC gauge prescrip- 
tion would give a framework in which a^js is fully determined by the target LC 
wavefunction, solved in the presence of an external field. 

8. Conclusions 

We have found that final state Coulomb rescattering in the target, within the 
coherence length 2v/Q'^ = 1/Mxb of the hard process, infiuences the iN — > i'X 
DIS cross section. In particular, diffractive scattering of the outgoing quark-pair 
on target spectators is a coherent effect which is not included in the light-front 
wave functions, even in light-cone gauge. Such effects modify the contributions 
of the individual target partons, implying that the DIS cross section is not fully 
given by the parton probabilities of the initial state. These coherent effects are 
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reminiscent of the LPM effect , which suppresses the bremsstrahlung of a high 
energy electron in matter due to Coulomb rescattering of the electron within the 
formation time of its radiated photons. 



Our analysis, when interpreted in frames with g"*" > 0, also supports the color 
dipole description of deep inelastic lepton scattering at small x. Even in the case 
of the aligned jet configurations, one can understand DIS as due to the coherent 
color gauge interactions of the incoming quark-pair state of the photon interacting 
first coherently and finally incoherently in the target. 

Our analysis in light-cone gauge resembles the "covariant parton model" of Land- 
shoff , Polkinghorne and Short p!8| , in the target rest frame. In this description 
of small X DIS, the virtual photon with positive first splits into the pair pi and 
P2- The aligned quark pi has no final state interactions. However, the antiquark 
line p2 can interact in the target with an effective energy s oc k'j_/x while staying 
close to its mass shell. Thus at small x and large s, the antiquark p2 line can 
first multiple scatter in the target via pomeron and Reggeon exchange, and then 
finally scatter inelastically or be annihilated. The DIS cross section can thus be 
written as an integral of the cr{qp — > X) cross section over the p2 virtuality. In 
this way, the diffractive scattering of the antiquark in the nucleus gives rise to 
the shadowing of the nuclear cross section a{qA — * X) 

Our results do not contradict the QCD factorization theorem for inclusive 
reactions in a general gauge. However, they show that the apparent equivalence 
between the DIS cross section and the target parton probabilities (|I|) suggested by 
the forward matrix element (§) in A'^ = gauge is incorrect. The A± components 
of the gauge field give leading twist contributions in LC gauge. 

Our investigation was triggered by the fact that the physically plausible and 
phenomenologically successful Glauber-Gribov description of DIS shadowing 
[T0[] implies that final state interactions influence the DIS cross section. The 
physics of shadowing is associated with final state diffractive scattering rather 
than with the (real) light-cone wave function of the target. There remains the 
possibility of incorporating shadowing in the target wave function by solving 
it under the specific boundary conditions implied by the Kovchegov LC gauge 
prescription |]T^. 

Our analysis is consistent with the standard Operator Product Expansion. Hence 
the usual sum rules of the parton distributions remain valid in spite of the rescat- 
tering (shadowing) physics. We have not estimated the quantitative importance 
of the rescattering effects on aois, but it is natural to expect that they are more 
prominent at small values of xb where the coherence length is long. In par- 
ticular, diffractive DIS is related to shadowing and is apparently generated by 
rescattering contributions. 
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APPENDIX 

A. One-loop calculation in = gauge. 

In this Appendix we present the calculation of the two-gluon exchange amplitude 
B (^5|) in light-cone n ■ A = A'^ = gauge of a scalar abelian theory. We shall 
take the target mass to be of the order of the transverse momenta, i.e., rather 
than ([19|) , we here consider the kinematic limit 

2z/ ~ > > ki^, p2±, k', m, M-:^ kf , k+ = Mxs+pt (53) 

and show that the expression for the amplitude remains the same. 

Leading contributions to the amplitude can come from diagrams Ba . . . B^. of Fig. 
7. The factors associated with the gluon propagators are approximated as 

{p + p'rd,,{k){2l + kf ~ [D{k + /) - Dil)] (54) 

where only the part of the propagator (11) contributes, and the function 



D{p) = D{p±) is defined in (0). Similarly, the factor from the four- leg scalar 
abelian vertex simplifies to 

{2p - k,rd,,{k,)d^'%k,){2p' + k,),. ^ (55) 

1 2 1 I 2 I 

where again the rf^^ components dominate. A factor 2ig'^ has been omitted for 
the time being. Direct use of the Feynman rules and of the kinematics (^) leads 
to: 

71 kik2 

[D{p,) - Dip2 - k2)][D{p2 - k2) - D{p2)] 



Ba + B, = -eg^QM I^I^Jik^)! 



-DiPi)[-P2kt + D{p2) - D{p2 - k2) + le] 
dki J,., f dk2 1 



B, + Ba = -eg'QMj^J^J{k^)J 



X 



[Dip^) - D{p2 - k,)][D{p2 - fci) - D{p2)] 



[-Plkt + D{pi) - D{p2 - ki) + ie][p2kt - D{p2 - k{) + le] 
dki j,,-^ [ dk2 1 ki± ■ k2 



- -es^QM I I ^ Jik^) I ^ ' ^ (56) 
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(a) 



(b) 




Figure 7: Diagrams that can give leading order contributions to the one-loop 
amplitude B in = gauge. 

where we use the shorthand notations 

d% I 1 



^^^1^ " + klJM -le^ + klJM - le ^^^^ 

In order to isolate the poles at kf = coming from the gluon propagators we 
view the integrands in (^) as rational functions of ^2", which we decompose in 
terms of simple elements. Also, since Pi is the largest scale we can approximate: 

1 1 / 1 1 \ 



ktl-plk^ + DiPi) - D{p2 - ki) + te] D{p^) - D{p2 - ki) \kt k+ - le ^ 

(5 

We also use 

1 1/1 l\ 



1 2 \ 1 



(59) 



2 
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to arrive at 

Ba + Bb 



-eg^'QMp, 



dki 

± J 2171 

1 



kt - \D{p2) - D{p2 - k2)]/p2 - le 

1 D{p2 - k^) 



Be + -Bd 



1 / _ D{p2-k2) \ _ 
kt\ D{p,) ) 



^2 



D{P2 



(60) 



dk:l 



D{p2) D{p2 - ki] 



1 



B. 



A^2 ^2 
,4, 



-eg'^QMp-^ 



dki 
± J 2i7r 



J(A;r 



dk2 —ki± ■ k2± 



1 1 



Using the relation 

2h± ■ k2± = D{p,) + D{p2) - D{p2 - k,) - D{p2 - k2) 
one easily checks that the terms oc 1/kf in (pO) give the contribution 



kt 



1 



dkt 

1 , D{p2 - ki) - D{p2 - k2) 



(61) 



(62) 



1 

^2" 



D(P2) /^(P2 - k2) D{pi 

1 2 



^(Pl)^(P2) 



D{p2) D{p2 - ki) D{pi 



1 ^ D(p2 - A;2) - D{p2 - ki) 



D{pi)D{p2) 



which vanishes by symmetry of /j_ and J{k^) under (A;^, /cf , kii_) (A;^, /c^, fc2±)- 
As a consequence, the sum of all diagrams is independent of the way one regu- 
larizes the spurious poles at kf = Noting that 



dki 

'2i7T 



Aki 



1 



the prescription independent result for B reads 

1 2 



B = -ieg'^QMp- 



+ 



(63) 



(64) 



D{p2) D{p2 - k2) D{pi) 

in agreement with the result (|35|) in Feynman gauge (and large M). 

As an individual diagram may contain pole terms ~ 1/kf, its value can depend 
on the prescription. As an illustration, we give the expressions of the different 
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diagrams using the three following prescriptions: 



{kt-^Vim + ^V^)V (PV) 



kf — ir], 
kf — ir]ie{k. 



(K) 

(ML) 



(65) 



namely the principal-value, Kovchegov^ and Mandelstam-Leibbrandt |]16 
prescriptions. The 'sign function' is denoted e{x) = 6(x) — 9(— x). With the PV 
prescription we have 





' 1 " 


/ dkt 









(66) 



V2 



and get 



B„ + Bj, 



-leg^QMp^ 



1 



1 



Be 

Using the K prescription we obtain 



-teg^QMp^ 




D{pi) D{p2 - k2) 

1 1 



D{p2) D{p2 - k,) 



(67) 



Ba + Bb = -2ieg^QMp^ 
Bc + Bd = 



1 



5, 



-leg^QMp-^ 



D{p^) D{p2 - k2] 
2ki± ■ k2i_ 



1 - 



D{P2 - k2) 

D{P2) 



(6^ 



D{pi)D{p2) 

The calculation with the ML prescription is a little more complicated. Defining 

h = j dk^J{k^)Q{-k2) 

I2 = f dk^J{k^)&{k^) (69) 



and using (|63|) we get after regularizing 

eg^QMp2 

ija + Db = 

X 



D{p^) D{p2 - k2 

\ D{P2 - k2) ( D{p2-k2) \ 



^Only the d ^ component of the gauge field propagator in Eq. (4) of contributes in our 
calculation. 
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Be + Bd 



B, 



TT 

eg^QMp^ 



1 



1 



TT 



D{p2) D{p2 - h] 



D{pi)D{p2] 



(70) 



k"^ 



Calculating explicitly Ji and I2 gives 

h = log 

h = log 
We can then use the relation 



IP 


-kl\ 


\kl- 


- h"^ 1 



+ Z7r[l + e(fci-A;2 
+ z7re(A;2_L 



ki) 



Ii{kuk2)+l2{k2,ki) = 2z7r 



(71) 



(72) 



to check that the sum of all diagrams evaluated with the ML prescription indeed 
reproduces the result (|5^). 

Instead of using (^), one can also directly use (^), after regularizing the kf = 
poles with a chosen prescription (for instance one of those given in (§^)), and 
perform the /c^ integral using Cauchy's theorem. The calculation is more in- 
volved, but reproduces all results presented above. See the comments at the end 
of Appendix B concerning this procedure. 

B. Two-loop calculation in = gauge. 

Here we evaluate the three-gluon exchange amplitude C (^7^ in = gauge and 
in the kinematic limit (|53D . The leading order diagrams Ca ■ ■ - Cg are displayed 
in Fig. 8. For each diagram, the 6 permutations of the vertices on the target line 
are taken into account. Since two permutations correspond to the same topology 
for diagrams Cd - ■ - Cg, there is a factor 1/2 for those diagrams. We will use the 
following shorthand notations: 



d'^ki_i_ <Pk2i_ 



(27r)2 (27r)2 kl^kl^k 



2 

■3± 



2tx 27r 



J 



K + kl^/M - ie A;f + k2 + {ki± + k2±)'^/M - ie 



dk-^ dk'2 
27r? 27ri 

perm(A;i, /csj^s) 
(73) 



(where J contains 6 terms arising from the 6 permutations mentioned above), 
and 

Ai = D{pi^ - kj±) ; A = D{pi^) ■ for ^ = 1, 2 and j = 1, 2, 3 (74) 
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Figure 8: Diagrams that can contribute to the two-loop amphtude C in = 
gauge. All permutations of the attachments to the target line are implied. 



where D is defined in (pTI). Using the kinematic limit (|53|) and approximations 
as in (|54D and (|55|) , the scalar abelian Feynman rules give: 

C - -2ea'0M I I J I ^ 

zeg J ^^^^^^ -D,{p^kt - + it){-p^kt - N; + le) 



X 



X 



1 



1 



= -2eg'QM j^j J 
Ce = -2eg^QM J J J 



-Pi {kt 


+ k+ 


) + Ni + N2 + ie){-p2k 




■ k2i_ 


1 




k^ 


^Di{-p2kt - N^ + te) 




■ h± 


1 




kt 


-Di{p2kt - Dn + te) 
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with 



-2e/QM 



u 









N3 h± ■ k2± 

ft,2 rt,3 

1 



-Piikt + kt) + + ie){-p2kt -N3 + le) 
9 ^DM f f T f ^1 • fc3± 1 



Ni = Di- ; iV2 = Dii - D23 ; N3 = D23 - 



(75) 
(76) 



Similarly to what was done in Appendix A for the one-loop calculation, one now 
considers all integrands in (^) as rational functions of kf {i = 1, 2, 3), which we 
decompose in simple elements, making first use of 



1 



1 1 

+ 



+ 



1 



th^ ryj'2 rv \ rv-^ rv^ '^3 '^2 3 



(77) 



The limit Pi ^ 00 must be taken after the decomposition in simple elements has 
been completed, otherwise some pinch singularities can arise. As there are in the 
two- loop case two independent '+' integration variables {kf + k2 + k^ = /c+), 
each integrand can be expressed as a sum of terms having one of the following 
forms: 

^ ^ ^ i^^j) (78) 



ktk+ ' kt{kf±ie) ' (fc+±ie)(fc+±z 



In (|78]) the poles at k^ = come from the gluon propagators in LC = 
gauge, whereas those at kf' = ±ie originate from the scalar quark propagators. 
The result of the full decomposition is 



Ca 



X 
+ 
+ 



-2eg^QMp^ J^J J 



\D23 

1 



1 \ Dn 



1 



DJ\D2 



- 1 



1 



D 



23 



X 



Dii DJ D2 

{Di - Dn){Dn - D23) 
DiD2{D2 - Dii) 

~2eg'QMp^ J^l J 
.D2 D23) V Dn 



1 



1 



k{ + ie k{ 
1 1 



1 

k2 



fcg — ie 



kj + ie kj 
1 ^ 1 



1 1 

fi,2 fi-S 



kt 



kt 



kt — ie 



kt — ie 



le 



le 
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+ 









" 1 1 " 




" 1 1 








^3 ^3 




kf — ie 



= -2eg'^QMp^ J^J J J 



X 



( ' 






-D23 ' 


\D23 


D2J 




-D23 



1 



kt + ie 



fci I ■ k 



2± 



D1D2 



1 

kt 

1 

T 



1 



1 

14~~ 



kt 



te 



+ 



1 - 



D 



23 



1 
kf 



kj + ze 



1 



^3 — ie 



+ 



-2eg'QMp^ J^l J 



ko\ ■ k 



3± 



DID2 



El 

D 



11 



1 1 

7T + T^ 



kj + ze A;]^ 



Ct 



D2 - Du 

-2eg'^QMp. 



1 

k^ 



1 



1 



/^23 



^2 
J 

1 



le 



1 



fcf + ie 



ki± ■ k 



2± 



Di - D23 
1 



ki — ie 
J 



1 



kt 

-k2± ■ ks 
D2Dn 



1 ^ 



1 



1 



1 



(79) 



Eq. (|7Up can be conveniently used to group together the poles at kf = 0, which 
appear in the two first forms of ([78|). For each of these forms, a lengthy calculation 
shows that the kf = poles add to a contribution which is identically zero, 
analogously to (|6^) for the one-loop calculation. On the way we use the identities 



- 2ki± ■ k2± 
-2k2± ■ h± 



Dis 



D12 
D12 



D23 
D21 



Di 
Di 



(80) 

D2 + Dn + D12 + Dr3 - D21 - D22 - D23 - D, = (81) 

and realize that in every factor [l/kf ± l/(^/ ± ie)] (i 7^ j) of ([79D, fc^ can be 
replaced by kf (and not the contrary^) by a change of variable. We also use the 
symmetry of and J under k^ ^ kj for i ^ j. 

Thus we have explicitly checked the complete prescription independence of our 
two- loop calculation. Only terms of the last form of (^) remain in (ffOl) . Using 



J J (k^ , k2 



12) 



^We do not allow the inverse change 1/kf V^/ to keep the possibility to deal with a 
regularized form of \/kf depending on ^ as is the case for the ML prescription, see (|65|). 
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as well as (0), ( |8lD and symmetry arguments, one shows that these terms add 
to 

3 3 



C 



1 



1 



D,, ' Z^n A J ^^^^ 

which exactly reproduces the result (^7|) obtained in Feynman gauge (for large 
M). 

After having shown the complete prescription independence of our calculation, we 
conclude this Appendix with some important remarks. We stress that (|75| ) and 
( |79[ ) are equivalent mathematical expressions for any of the diagrams Ca- ■ - Cg. 
To evaluate a given diagram, one needs to regularize the kf = poles, but this 
can be done starting either from (^) or from (^), and the same results must 
follow. We have checked this for all diagrams using the PV and K prescriptions. 
We thus see no problems in applying the PV prescription to two-loop diagrams. 
Using the PV prescription on (^) is straightforward, but applying it to ([f5|) 
requires some comments. Regularizing l/{kik2k^) yields 



1 



n^l 



npv 




lim lim lim 

m^O r)2^0 r]i^O 



1 




1 




1 




[kt\ 






r)2 


_k^ _ 


m 



^4) 



where [\-/kf]r,^ is given in (|^). Thus the poles at kf = must be regularized with 
distinct small finite parameters 77,. Then the kf integrals are performed using 
Cauchy's theorem, and only in the end the limits ?7i — >■ 0, 772 0, 773 — >■ are 
taken separately (in arbitrary order). We found this procedure to be well-defined 
and to give results consistent with those directly obtained from (ffT 



Finally, as in the one-loop case, it is remarkable that the K prescription makes 
all two-loop diagrams where the fast quark rescatters vanish, i.e., only Ca, Cd 
and Ce contribute to the amplitude C. 



C. Absence of struck quark rescattering in = gauge. 

In this Appendix we show that in ^4+ = gauge, independently of the prescription 
used to regularize the spurious kf = poles, rescatterings of the struck quark pi 
cancel in the cross section, i.e. after summing over cuts in the forward Compton 
amplitude. This is done by proving that the full contribution to the cross section 
(use Eqs. (P, m, (H) 



\B\ 



2AC 



--{eg'QMp2 



d^r^ d^R^ V{m\\r^fW{r^, R^ 



^5) 



is given by attachments to p2 only. 



We need to know the partial amplitudes A2, B2, C2 contributing to A, B, C 
where only attachments to p2 are kept. For the Born amplitude A, only the 
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diagram of Fig. 4a contributes in ^4+ = gauge. Thus the partial amphtude A2 
from attachments to p2 is actually the full amplitude A given in (p2D , 



A, = A 



2eg^MQp2 



1 1 



(86) 



The partial one-loop amplitude B2 is given by the sum of the diagrams Ba, Bj, 
and Be of Fig. 7. This sum is prescription dependent. We use the notation (see 
Eq. m) 



f dH 


' 1 " 


J in 





(87) 



J r) 



giving Ipv = and Ik = ^- After regularizing (^) and using (p3| ) and symmetry 
arguments we findQ, 



Bo = Bn -\- Bh + Bp 



Di D 



22 



J 1_ 

^'[02 D22 



The partial two- loop amplitude C2 is given by the diagrams Ca, Cd and Cg of 
Fig. 8. Regularizing ([79| ) and using (p^ ) we get 



C. ^ -ie/<3Mpj/^{(J--J-)(l 



D 



11 



23 



1 a. + 1) 



+ 



-eg^QMp-^ 



1 A Dn-%,^2_T^ C^WW^n^iMrr 



/^2 

ki± ■ k2± 
D1D2 

k2± ■ k3i_ 
D1D2 



2/^(/^ - 1) + ( 1 



DiD2{D2-Dn 
D2 



1 



D2Z 

2Ir,{-Ir, - 1 



Du-Di 
D2 - Dn 



{Ir, - 1) 



(89) 



Using (^0]), (^Tj) and symmetry arguments we get after some algebra 
C2 = Ca + Cd + Ce = -^eg''QMp, 



2 1 

+ 



3Di 6D2 2D 



1 1 

+ 



23 



D 



11 



+ 4 



D 



1 1 

+ 



23 



1 



1 



2D2 2D 



23 



(90) 



It is an easy exercise to express the partial amplitudes A2, B2, C2 in transverse 
coordinate space, as was done for the full amplitudes A, B, C in section 5 (see Eq. 

''In order to use (^3|) in (|60|), we need to consider a regularized form of l/k^ independent 
of k~ {i.e. we exclude for simplicity the ML prescription in this Appendix). Eq. (|8^) is valid 
for any such prescription. 



32 



). Since the partial amphtudes are not infrared finite^, we introduce a small 
photon mass A in the exchanged photon propagators, i.e. l/kfj_ — ^ + '^^) 

in the definition of { ^7\j or (73). Then 



A2 
B2 



2eg^QMp^VW 
-i eg^QMp^VW 
1 



I,Ko{XR^)-Ko{X\R^ + r^\) 



4 



eg^'QMp^VW 



2n 

( l,K^{\Ri)-K,{\\R^ + r^\ 
3 27r 



(91) 



where V and W stand for V^(m||r_|_) and Vr(rj_,i?^). The contribution from 
attachments to p2 to the cross section reads 



d^r^d^R^ \B2? + 2A2C2 



(eg^QMp^)^ d^f^d^R^V^W^ (92) 



This is infrared finite, prescription independent, and identical to the full result 
(pSl). Hence rescatterings of the struck quark pi cancel in the cross section in 
y4+ = gauge. 



From (85) and (92) we see that in = gauge and in coordinate space, the 
contribution from attachments to p2 equals the full contribution |-Bp + 2 AC even 
at the integrand level, i.e. before integrating over and R±. We thus have, in 
coordinate space, 

\M\^=\M2\^ (93) 

where M = A + B + C + ... is given in (39) and M2 = A2 + B2 + C2 + ■ ■ ■ 
corresponds to the partial amplitude where only attachments to p2 are kept. 
Thus in A'^ = gauge 



M{f±,R^)=e''^M2if±,Ri.) 



(94) 



The full amplitude Ai is obtained from 7Vl2 by inserting any number of rescat- 
terings of the quark pi. Eq. (R3) reads 



A+B+C+ 



e"^ {A2 + B2 + C2 + .. .) 



(95) 



By expanding the l.h.s. and r.h.s. of ( P5| ) up to order g\ one realizes that (j) must 
be at least of order (7^, 



hg'' + ... 



(96) 



Note however that with the K prescription (/^ = 1) the partial amphtudes ( |91| ) equal the 
full ones, as already mentioned at the end of Appendix B. Thus the partial amplitudes are finite 
when A with this particular prescription. 
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Identification of tlie terms of order and in the two sides of ( pSj) leads to 



4tt 







(97) 



or 



g' IlL^ KoiXR^) + Oig' 



(9^ 



Although not proven here, the O {g^) terms in are expected to vanish because 
adding one rescattering of pi can only bring a power g"^ (see also the following 
discussion). Thus we get 



M{f±, R±) = exp 



t^{Ir,-l)Ko{XRA 



M2(f±,RA 



(99) 



As expected, since M. is infrared safe and prescription independent, all the de- 
pendence on A and of A^2 is contained in the phase. Note also that with the 
K prescription. If, = I and M. = M.2- 

Eq. (^) can also be understood as follows. In momentum space, if we call /i the 
Lorentz index associated to the coupling of pi, we know that the amplitude is 
dominated by the n^k^ /kf term of the exchanged gluon propagator in = 
gauge, with =± and u = Together with the scalar quark propagator 



Aj oc {pi — h) — m + ^ —pi k^ + ai + is 
where = D{pi) — D{pi — ki), the factor l/k^' yields 



(100) 



dkt 



27ii (— Pi k:^' + flj + ie)k^ 



kt 



+ 



1 



1 /. 



J -q 



-kj + ie 



(101) 

where some simplification similar to (58) was made. The scalar coupling brings 
a factor g'^ki± ■ {2pi — ki)± = g'^ai which compensates the prefactor in the r.h.s of 
( |101| ). We are left with the l/(/cy_ + A^) factor from the gluon propagator, which 
after Fourier transform gives Kq{XR^) / {27i) . This builds the complete phase in 
191). 
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